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We consider coupled nonlinear Schrodinger equations (CNLSE) which govern the propagation 
of nonlinear waves in bimodal optical fibers. The nonlinear transform of a dual-frequency signal is 
used to generate an ultra-short-pulse train. To predict the energy and width of pulses in the train, 
we derive three new types of travelling periodic-wave solutions, using the Hirota bilinear method. 
[ We also show that all the previously reported periodic wave solutions of CNLSE can be derived in 

a systematic way, using the Hirota method. 
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PLh ■ I. INTRODUCTION 

d [ 

Steady increase of bit rates provided by fiber-optic telecommunication networks is necessary to support the rapid 
growth of the data traffic. The development of high-repetition-rate generators of high-quality ultra-short pulses 
plays a crucial role for the increase of transmission capacity in fiber-optic networks, as well as for optical-signal 
I , processing [1]. Since available opto-electronic direct-modulation techniques cannot provide for the generation of 
£>. ' pulses at a repetition rate higher than 40 GHz, nonlinear transformation of a dual-frequency signal in an optical fiber 
, is an attractive method for applications to ultrahigh-speed telecommunications. Indeed, the repetition rate of the 
t-H ■ generated pulses can be tuned by adjusting frequencies of the two input waves. A nonlinear technique based on the 
modulational instability of a pump wave induced by a small periodic signal was first suggested by Hasegawa [2] and 
realized experimentally by Tai et al. [3] . Another nonlinear method is based on reshaping of a beat signal which results 
from superposition of two pump waves of equal power. Nonlinear reshaping of the beat signal into well-separated 
short pulses was previously demonstrated by several methods, using dispersion decreasing fibers (DDF) [4], intra-pulse 
stimulated Raman scattering in DDF [5] , or switching in a nonlinear fiber-loop mirror [6] . Other techniques based on 
step-like and comb-like dispersion-profiled fibers, composed of segments of conventional uniform fibers with different 
dispersions, have also been proposed for the generation of picosecond and sub-picosecond pulse trains with high 
repetition rates [7-11]. However, the common feature of these techniques, employed for the reshaping of a beat signal, 
is that they rely on careful control of dispersion distribution in complex fiber systems. An alternative approach, based 
on multiple four-wave mixing (FWM) in an optical fiber with constant anomalous dispersion, has been demonstrated 
theoretically for the generation of compressed pulse trains from a dual-frequency pump field [12]. In particular, the 
optimum fiber length (propagation distance at which the pulse width attains a minimum) for given values of input 
beat-signal parameters, has been accurately predicted by calculating the propagation distance at which the maximum 
energy conversion from the pumps into the nearest FWM-sidebands occurs [12]. Experimental demonstration of the 
generation of a transform-limited, pedestal-free pulse train with high repetition rate by means of this technique of the 
multi- wave- mixing temporal compression in the standard non-zero dispersion-shifted fiber was reported recently [13]. 

For this kind of ultra-short pulse-train generation from the beating of two-mode signals, one must investigate 
periodic-wave solutions of nonlinear equations governing the fiber system. The single-mode wave propagation in 
fibers is governed by the well-known nonlinear Schrodinger (NLS) equation [14,15], 
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where tp is the slowly varying envelope of the axial electrical field and 2/3 is the group- velocity-dispersion (GVD) coef- 
ficient, while z and t are spatial and temporal variables. The normal and anomalous GVD corresponds, respectively, 
to (3 > and (3 < 0. 

Single-mode fibers actually support co-propagation of two different modes with orthogonal polarizations, which, in 
the case of anomalous GVD and without birefringence, is described by coupled NLS equations (CNLSE) [16-19]. 

iipu - Witt + + (IV-il 2 + ct|V 2 | 2 ) i>i = o, 

iihz - Witt + + (|^ 2 | 2 + <j\^i\ 2 ) fa = 0, (1) 

where a > is the relative cross-phase modulation (XPM) coefficient, and the linear coupling coefficient n accounts for 
possible twist of the fiber in the case when ipi and tp 2 represent orthogonal linear polarizations, or elliptic deformation 
of the fiber's cross section if ipi and -02 correspond to circular polarizations. In fact, the linear coupling will play a 
crucially important role in the present work. 

Complete integrability of the system (1) by means of the inverse scattering transform (1ST) in the case a = 1, k = 
was shown by Manakov [16]. In a completely integrable system, one can use 1ST to construct not only solitons, but 
periodic solutions too. However, a ^ 1 in real optical fibers, hence Eqs. (1) arc not integrable. In particular, a = 2 
for two modes with circular polarizations, as well as for a different situation, when two signals carried by different 
wavelengths co-propagate in the fiber. For two linearly polarized modes, one has a = 2/3, and in the more general 
case of clliptically polarized modes, the XPM coefficient may take any value from 2/3 < a < 2. All these physical 
situations come under the umbrella of the non-integrable CNLSE system. 

The objective of the present paper is to find new periodic-wave solutions of the CNLSE system (1) for both 
anomalous-GVD and normal-GVD regimes. Solutions will be obtained by dint of the Hirota method (HM) [20] 
through the Jacobi's elliptic functions. Links to previously found periodic-wave solutions of the same system will be 
established too. 

The rest of the paper is organized as follows. In section 2, we produce three new types of periodic-wave solutions 
for the anomalous-GVD case, and two new solution types for the normal-GVD case. In section 3, we briefly consider 
the soliton (long- wave) limit of these solutions, concluding that they amount to anti-dark and dark (generally, gray) 
solitons in the anomalous- and normal-GVD cases, respectively. Conclusions are formulated in section 4, and some 
technicalities are collected in two Appendices. 



II. PERIODIC WAVES IN THE CNLSE SYSTEM 



A. The anomalous-dispersion regime 

In this subsection we consider the case (3 < in Eqs. (1), adopting the normalization (3 = — 1. Some periodic- wave 
solutions of nonintegrable CNLSE systems have been found earlier for this case in Refs. [21-23]. Here we shall obtain 
three new kinds of periodic solutions by means of HM, in terms of the theta and Jacobi elliptic functions. We shall 
also show how previously known periodic solutions of the system (1) can be readily found by means of HM. 

We start with the transformation 

il> 1 = ^(a + ib), 4> 2 = ^{a- ib) , (2) 

which casts Eqs. (1) with (3 = — 1 in the form 

ia z + a tt + kg+ i(|a| 2 + |6| 2 )a + i(cr - l)(a 2 + 6 2 )a* = 0, 

ib z + b„ -Kb+ l -(\a\ 2 + \b\ 2 )b + 1((7 - l)(a 2 + 6 2 )6* = 0, (3) 

where the asterisk stands for the complex conjugation. Stationary solutions to these equations are sought for in the 
form 

a = exp(— iflz), b = rj(t) exp(— iSlz), 

where the functions £ and rj are assumed real. Coupled equations for these functions are tantamount to the stationary 
version of the Manakov system, 
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Z xx + (n + k)£ + -(1 + a)(£ 2 + ?7 2 K = o, 

77^ + (fi - K)r? + ^(1 + ct)(£ 2 + ry 2 )^ = 0. 



(4) 



Despite the asymmetry of equations (4), accounted for by the terms proportional to ±k, this system is an integrable 
one (it has a second dynamical invariant, besides the Hamiltonian) [17]. 
We now employ the Hirota bilinear operator: 
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dx' J 




dt> j 



g{x,t)f{x',t' 



(5) 



to rewrite Eqs. (4) in the trilinear form [24]: 



x—x'.t—t' 



t_9 _ G 



f(D 2 x + n + K )g-f + g 
f(D 2 x + n- K )G-f + G 



-Dlf.f+-(l + a)(g 2 + G 2 ) 
-Dlf-f+^l + a^ + G 2 ) 



0, 
0. 



(6) 



Periodic waves are obtained by choosing g, G, and / as products of theta functions. The Hirota bilinear forms of 
the theta functions can be simplified by using known properties of these functions. Technical details are omitted in 
the main text, as the method is very similar to that presented in Refs. [24,25]; a brief sketch of the derivation of the 
solutions displayed below is given in Appendix A. The results are presented in terms of Jacobi elliptic functions, rather 
than theta-functions, as the Jacobi functions are more straightforward to handle by means of computer-assisted algebra 
software. As an independent check we verify that all the solutions given below satisfy Eqs. (1) by direct differentiation 
with the software package Mathematica. 

An important point is the choice of the theta function for /. The functions $i(x) and $2(2:) have real zeros, therefore 
choosing them will give rise to singular solutions. The functions #3(0;) and $4(2) are related by a phase shift of tt/2. 
Previously known solutions [21-23] employ the choice / = "&\{x). In this work, we start with / = #3 (a;) #4 (a;), hence 
the results are different from the previously known ones. Furthermore, the advantage of HM is that it can be readily 
generalized for CNLSE with more than two components. 

By means of this technique (see details in Appendix A), we have obtained three new sets of solutions for the system 
(3). The first set has the form 



3 <A 



^2 = 2^1 



(1 - k 2 ) lfi dn(rt) 



dn(rt) (1 - k 2 ) 1 ^ 
(1 - fc 2 ) 1 /* dn(rt) 



+ ikB lS n(rt)cn(rt) \ c _ iQ 



dn(rt) 
ikBisn(rt)cn(rt) 
dn(rt) 



dn(rt) (l-fc 2 ) 1 /4_ 

where k is the modulus of the elliptic functions. The amplitude parameters in these expressions are related by 

fc 2 ,4 2 8r 2 fc 2 



B\ = 



vT^fc 2 l + cr 



the frequency is 



n = rVl - k 2 + ^r 2 (2 - k 2 ) - ^A\{\ + a) 



Vl-fc 2 



1 



a/I - fc 2 



and the modulus k is related to the linear coupling constant n: 



K = r 2 (Vl — k 2 - 1 + ifc 1 



(7) 



(8) 



(9) 



(10) 



Thus, the solution (7) contains two free parameters, as five constants Ai,Bi, k, r, and are subject to three constraints 
(8), (9), and (10). Each of two other solutions displayed below also contains five constants on which three conditions 
are imposed. An example of the periodic solution (7) is shown in Fig. 1 (a). 
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The second solution is 



(l-fc 2 ) 1 /* 



dn(rt) 



^2 = g 1 



dn(rt) (l-fc 2 ) 1 / 4 
(1-fc 2 ) 1 / 4 ^ dn(rt) 



dn(rt) (l-fc 2 ) 1 / 4 



zfc.Bisn(rf)cn(rt) ■ e _ ir , ; 

dn(rt) 
ikB lS n(rt)cn(rt) \ _ in 

dn(rt) 



(11) 



The relation between the amplitude parameters in this solution is the same as given by Eq. (8), but relations between 
other parameters are different, 

SI = -rVl — fc 2 + ^r 2 (2 - fc 2 ) - ]a 2 (1 + a) y/l - fc 2 + - =L + 



2 v '4 

n = -r 2 fc 2 + 1 - ^k 2 ^j 

Figure 1 (b) shows an example of this solution. 



yT^fc 2 



(12) 




FIG. 1. (a) The intensity \ipi\ 2 versus t for the solution (7) with the amplitude relation (8) in the case Bi = 1. (b) The 
solution (11) in the case B\ — 1. (c) The solution (13) in the case A\ — A-2- In all the three cases shown, and also below in 
Fig. 2, r = 1 and fc = 0.7. Note the difference in the vertical scale between the three solutions. 



The third explicit solution is 



V>i = 2 i A i 



*h = 2 i M 



(l-fc 2 ) 1 / 4 dn(rf) 



dn(rt) (l-fc 2 ) 1 / 4 
(l-fc 2 ) 1 / 4 dn(rf) 



dn(rt) (1 - fc 2 ) 1 / 4 



+ iA 2 
- iA 2 



(l-fc 2 ) 1 / 4 dn(rt) 



dn(rt) (l-fc 2 ) 1 / 4 
(1-fc 2 ) 1 / 4 ^ dn(rt) 



dn(rt) (l-fc 2 ) 1 / 4 



-iQz 



iilz 



(13) 
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This time, the amplitude parameters are related by the equation [cf. the relation (8) for the two previous solutions] 



(2 , a2 8r 2 VT^W 



A( + Ai = ^- , (14) 

1 + cr 



and other constants obey relations 



Q = -r 2 (2 - fc 2 ) - 4rVl - fc 2 + (1 + a)A 2 , 



k = 2rVl - k 2 . (15) 

Figure 1 (c) shows an example of this periodic solution. 

All the above solutions exist only in the presence of the linear coupling: as it follows from Eqs. (10), (12), and (15), 
no nontrivial solution can be found if k = 0. On the other hand, the solutions exist at any value of a, the Manakov's 
case, cr=l, having nothing peculiar. 

As it was mentioned above, the same HM-based procedure makes it possible to reproduce all the other periodic 
solutions of the CNLSE system (1) which were reported earlier in Refs. [21-23]. Those solutions can be obtained on 
the basis of different combinations of the theta functions, see details in Appendix B. 



B. The normal-dispersion case 

Proceeding to Eqs. (1) with > 0, we adopt the normalization = 1, thus dealing with the equations 

tyiz ~ i>itt + >#2 + (IV-il 2 + <?\H 2 ) V>i = 0, 

iip2 Z - ihtt + kiPi + (\H 2 + ^IV-il 2 ) V>2 = 0. (16) 

As well as in the case of the anomalous GVD considered above, we can repeat the same steps, involving the trans- 
formation (2) and HM, to generate periodic-wave solutions for the normal-GVD regime. We have thus obtained two 
new sets of solutions, using different identities for theta- functions. 

The form of the first solution is the same as given by Eqs. (7), but the amplitude parameters are related by 



The relations (9) and (10) for the solution (7) are replaced by 

n = -\r 2 {2 fc 2 ) - r 2 + o)A\ [^k 2 + 2) , (18) 



..2 



1 - h 2 - yA~k^j . 



(19) 



The form of the second solution is the same as given by Eq. (11), and the amplitude parameters are related by 
(17). However, the other parameters are related by 

O = -^r 2 (2 - fc 2 ) + rVl-fc 2 - \{1 + <J)A\ (^1-k 2 + -j=L= + 2) , (20) 



K = r 2 (V 1 - k 2 + 1 - ifc ; 



(21) 



[cf. Eq. (12)]. Figure 2 shows an example of the latter periodic solutions. 
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FIG. 2. The intensity \ipi\ 2 versus t for the solution (7) of the dark-soliton type, subject to the amplitude relation (8), in the 
case B\ = 1. 

There is no counterpart for the third anomalous dispersion regime solution (13) in the normal-GVD regime. 



III. THE SOLITON LIMIT 



A natural question is the soliton (long-wave) limit of the periodic solutions derived above for both the anomalous- 
and normal-GVD regimes, which corresponds to k — > 1. Note that many expressions in the above solutions seem 
singular in this limit [for instance, Eq. (14)]. However, accurate analysis shows that a nonsingular limit exists. 

In the case of anomalous GVD, both solutions (7) and (11) give rise to a single long-wave limit, 



1p2 



8r 2 

f3f + sech (rt) + zfttanh (rt) 



l + <7 



0i 



8r 2 
T+a- 



sech (rt) — i/Jitanh (rt) 



exp 



exp 



where is an arbitrary real constant, and 



K=--r . 



(22) 



(23) 



The solution (22) corresponds to a state in the form of a bright soliton placed on top of a continuous- wave background, 
as it is seen from the expression for the intensity: 



l^ll 2 = \lp2 



8r 2 

01 + sech 2 (rt) 

1 + a 



Solutions of this type are often called anti-dark solitons. 

The long-wave limit of the solution (13) is only possible, as it follows from Eq. (15), if the linear coupling n vanishes 
as \J\ — k 2 simultaneously with (l — fc 2 ) [on the contrary to the above soliton- limit solution, see Eq. (23)]. In the 
corresponding limit form of the solution, both ipi and ^2 are proportional to sech (rt), which implies an obvious bright 
vector soliton. 

In the case of normal GVD, the common long-wave limit of the solutions given by Eqs. (7), (17), (18), (19) and 
Eqs. (11), (17), (20), (21) is 



"01 = 



-fosech (rt) +i\ 0% 



8r 2 

l + CT 



tanh (rt) 



8r 2 



-/3 2 sech (rt) - i\ 0% + tanh (rt) 

V 1 + cr 

where /?2 is an arbitrary real constant, and 



cxp(-m { s 2 2z), 



exp(-i^ 1 ) z), 



\r*-\(l + v)Pl 



K=-r . 
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As expected, this is a dark soliton with the intensity distribution 



m 2 = -M + 



8r 2 



tanh 2 (ri) . 



l + a 



In the general case (/?2 ^ 0), this dark soliton is of the gray type, i.e., with a nonzero minimum value of the intensity 
at the central point, t = 0. 



The periodic solutions obtained in this work for the nonlinear fiber-optic model with both anomalous and normal 
GVD can be helpful in predicting the energy and width of ultra-short soliton trains using the bi-modal pulse prop- 
agation in nonlinear optical fibers (the underlying CNLSE model is definitely relevant for the pulse trains with the 
temporal width of the pulses up to <~ 1 ps, provided that the loss is compensated by appropriate gain). Using these 
analytical solutions, and knowing fiber parameters, such as GVD, SPM and XPM coefficients, one can calculate the 
exact shape, width and power of the periodic solution. These results may also help to perform numerical simulations 
of the beatings of dual-frequency signals. 

The above consideration was not dealing with stability of the pulse trains. In the case of the anomalous GVD, 
the trains may be unstable; however, in many cases one may expect that their instability is quite weak (see, e.g., 
Ref. [26]), which makes it possible to consider practical applications of these waves, especially in the case of limited 
propagation distance, which will not give the instability enough room to develop. On the other hand, the pulse 
train may be completely stable in the case of normal GVD. Moreover, a stable wavetrain with the carrier frequency 
chosen in the normal-dispersion spectral region can be used to stabilize a co-propagating periodic wave launched in 
the anomalous-GVD region (on the other side of the zero-dispersion point), provided that group velocities of both 
waves coincide, as it was proposed in Ref. [27]. In fact, it is a challenging issue (which will be considered elsewhere) 
to find exact solutions for periodic waves in the bimodal system, with the XPM coefficient a — 2 and without linear 
coupling (k = 0), in the case of opposite signs of the dispersion in the two modes. 

For completeness, we have also presented the long- wave limit of the newly obtained periodic wave solutions, in both 
the anomalous- and normal-GVD cases. In the former case, it turns out to be an anti-dark soliton, while in the latter 
situation the long-wave limit amounts to a dark (generally, gray) soliton. 

To summarize, we have considered the CNLSE model (1), which governs the propagation of electromagnetic waves 
in nonlinear bimodal optical fibers. The model includes both XPM (nonlinear) and linear couplings, the latter 
one accounting for a twist or elliptic deformation of the fiber, in the cases of the linear and circular polarizations, 
respectively. Using properties of elliptic theta functions and the Hirota bilinear method, we have found three new 
types of periodic solutions in the case of anomalous GVD, and two new types in the case of normal GVD. The new 
solutions require the presence of the linear coupling. We have also shown that all the previously reported periodic 
solutions of the CNLSE model can also be obtained by means of the Hirota bilinear method. Thus, we conclude that 
this method furnishes a powerful tool to derive different forms of periodic-wave solutions in terms of elliptic functions. 
Basically the periodic solutions derived using the elliptic functions are in the form of raised sine or cosine functions, 
as seen from Figs. 1 and 2. In most of the return-to-zero (RZ) coding the commonly used optical modulators utilize 
Mach-Zender interferomctric technique. This technique yields raised cosine form of pulses for RZ bits coding. Hence 
we believe that the periodic solutions derived using Jacobi elliptic functions are of the same type as those generated by 
the Mach-Zender interferometer modulators. New analytical results obtained in this work may provide for important 
information in determining properties of pulse trains that are generated by dual-frequency signals in nonlinear optical 
fibers. 
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APPENDIX A: THE THETA-FUNCTIONS IDENTITIES 



The theta functions 9 n (x), n = 1,2,3,4, are denned by as follows below, together with auxiliary parameters q and r: 

oo 

0i (a;) = 2 J2(- 1 ) n Q {n+1/2)2 sin [(2n + !)&] 



n=0 

oo 

= -^cxp 

— OO 

oo 

2 (*)=2]T ( 

71=0 



:t / 7 i m + - 



2i[m + 



(n+l/2y 



cos [(2n+ l)x] 



J^exp 



7HT m + 



+ 2i[m + 



Inx) 



h{x) = l + 2^q n cos(2r 

n=l 

oo 

= ^ exp(7rirm 2 + 2imx), 



$ 4 {x) = l + 2^(-l)"g" 2 cos(2nx) 



71=1 



oo 

= ^exp[ 



7rirm + 2im 



q = cxp(7rir) 



(Al) 



Note that, in the present case, r is imaginary, hence the constant q defined in Eq. (Al) belongs to the interval 
< q < 1. It is related to the modulus k of the corresponding Jacobi functions as q = cxp [— ttK — k 2 ) /K(kjj , 
where K is the complete elliptic integral of the first kind. Zeros of the functions 6i, 62, #3, and 84 are at the points 
x = Mtt + Nttt, (M + 1/2)tt + Nttt, (M + l/2)n + (N + 1/2)ttt, and Mn + (N + 1/2)ttt, respectively, where M and 
N are integers. 

There are many identities involving products of theta functions which are relevant to the Hirota bilinear method. 
An example is 

+ y)Mx - y)02(o) - vl{x)#l{v) + 4(^i(y)- (A2) 

A proof can be found in Ref. [25]. On differentiating Eq. (A2) with respect to y twice and setting y = 0, one obtains 

2^(0)^ (*) 



DiMx) ■ 4 (aO = 2^(O)t?2(O)^(a0 



02(0) 



(A3) 



where the Hirota derivative [see Eq. (5)] is 

D 2 x g ■ f = fg X x - 2.fx9x + gfxx- 

From these results, one can construct equations of increasing complexity by using Hirota identities like 

D 2 x ab ■ cd = bd{D 2 x a ■ c) + ac(D 2 x b ■ d) + 2{D x a ■ c)(D x b ■ d). 

The novel choice in the present paper is to set / = 03(cc)04(a:) as opposed to the selection of / = 0|(x) used implicitly 
in earlier works (see Appendix B). The relevant identity is then 

Dlf ■ I = D 2 J 3 (x)Mx) ■ M*)Mx) 

= 2mm [*$(*) + + gg^MM , 

The trick is then to select the amplitude parameters such that only the term 0|(a;)0|(a;) remains in the square bracket 
of Eq. (6). The remaining terms of the trilincar equations are then satisfied by adjusting the frequency parameter f2. 
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APPENDIX B: COMPARISON WITH PREVIOUSLY KNOWN RESULTS 



Previously known periodic solutions of the CNLSE model (1) can be recovered by the present approach via proper 
choice of theta functions in Eqs. (6). For instance, selecting 

/ = 2 4 (at), 

24a 2 tf| (0)02(0) 



g = A 6 1 (at)6 2 {at), A z 



(1 + <7)0§(O) 



r r o( +\a< A r 2 24^(0)^(0) 
reproduces a solution reported in Ref. [22], which takes the following form in the present notation: 



a \ -j— ^ — rk 2 sn(rt) cn(rt) cxp(— iflz), (Bl) 



b = ^/y^— rkcn(rt)dn(rt)exp(-iQz), (B2) 



0=^ 2 (l-2fc 2 ), n=\r\ 

a is the wavenumber in the theta function representation. Since arguments of theta and elliptic functions differ by 
a scalar (Eq. (B3) below), the wavenumbers in the two formulations are related by r = a#§(0). The Jacobi elliptic 
functions appearing in the expressions (Bl) and (B2) are related to the theta functions as follows: 



dnfrt) - rt = -E- (B3) 

Note a difference in the scale factors in the arguments of the two types of elliptic functions. 
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